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ABSTRACT 

We present a hard X-ray spectrum of unprecedented quality of the Galactic supernova remnant 
W49B obtained with the Suzaku satellite. The spectrum exhibits an unusual structure consisting 
of a saw-edged bump above 8 keV. This bump cannot be explained by any combination of high- 
temperature plasmas in ionization equilibrium. We firmly conclude that this bump is caused by the 
strong radiative recombination continuum (RRC) of iron, detected for the first time in a supernova 
remnant. The electron temperature derived from the bremsstrahlung continuum shape and the slope 
of the RRC is ~1.5 keV. On the other hand, the ionization temperature derived from the observed 
intensity ratios between the RRC and Ka lines of iron is ~2.7 keV. These results indicate that the 
plasma is in a highly overionized state. Volume emission measures independently determined from 
the fluxes of the thermal and RRC components are consistent with each other, suggesting the same 
origin of these components. 

Subject headings: ISM: individual (W49B) — supernova remnants — radiation mechanisms: thermal 
— X-rays: ISM 



1. INTRODUCTION 

W49B (G43.3-0.2) is a Galactic supernova remnant 
(SNR) with strong X-ray line emissions from highly ion- 
ized atoms. It exhibits centrally filled X-rays inside a 
bright radio shell with a radius of 100 arcsec (Pye et al. 
1984). The distance to W49B remains uncertain. It was 
estimated to lie at a distance of ^8 kpc (Radhakrish- 
nan et al. 1972; Moffett & Reynolds 1994). Using the 
spectral distribution of HI absorption, however, Brogan 
& Troland (2001) found no clear evidence that W49B is 
closer to the sun than W49A, which is located at a dis- 
tance of ~11.4 kpc (Gwinn et al. 1992). This may indi- 
cate a possible association of W49B with the star- forming 
region W49A. The near-infrared narrowband images in- 
dicate a barrel-shaped structure with coaxial rings, which 
is suggestive of bipolar wind structures surrounding mas- 
sive stars (Keohane et al. 2007). They also showed an 
X-ray image from Chandra, which has a jet-like struc- 
ture along the axis of the barrel. They interpreted these 
findings as evidence that W49B had exploded inside a 
wind-blown bubble in a dense molecular cloud. 

Using ASCA data, Hwang et al. (2000) showed that 
broadband modeling of the remnant's spectrum required 
two thermal components (0.2 keV and 2 keV) and signif- 
icant overabundances of Si, S, Ar, Ca, Fe, and Ni. They 
confirmed that most of the X-ray emitting plasma was 
nearly in coUisional ionization equilibrium (GIF). They 
also found evidence for Cr and Mn Ka emission. 

Kawasaki et al. (2005) claimed the presence of "overi- 
onized" plasma in W49B through the analysis of ASCA 
2.75-6.0 keV spectrum. They measured intensity ratios 
of the H-like Ka (hereafter Lya) to He-like Ka (Hea) 
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lines of Ar and Ca to obtain the ionization temperature 
(kTz), and found that kTz (^^2.5 keV) is significantly 
higher than the electron temperature (fcTc) determined 
from the bremsstrahlung continuum shape (^1.8 keV). 
Miceli et al. (2006) adopted the same analysis procedure 
for the XMM-Newton spectrum of the central region but 
found no evidence for the overionized state. 

In this letter, we report the firm evidence for overion- 
ized plasma in W49B using data from the X-ray Imaging 
Spectrometers (XIS: Koyama et al. 2007) onboard the 
Suzaku satellite (Mitsuda et al. 2007). 

2. OBSERVATION AND DATA REDUCTION 

We observed W49B with Suzaku on 2009 March 29 and 
31 (observation IDs 503084010 and 504035010, respec- 
tively). Because the pointing positions and rotation an- 
gles do not differ significantly, we merged these two data 
sets. The XIS consists of four X-ray CCD camera sys- 
tems placed on the focal planes of four X-ray telescopes 
(Serlemitsos et al. 2007). Three of them were operated 
in the normal full-frame clocking mode with a spaced- 
row charge injection technique (Uchiyama et al. 2009) 
during our observations. Two are front-illuminated (FI) 
CCDs and the other is a back-illuminated (BI) CCD. 
Data were cleaned using processing version 2.2.11.24. 
We used HEASOFT version 6.5.1 for data reduction and 
XSPEC version 11.3.2 for spectral analysis. After screen- 
ing with the standard criteria, ^ the net integration time 
was ^^113 ksec. 

3. ANALYSIS 

3.1. Construction of the Spectrum 

Figure [1] shows the vignetting-corrected XIS image in 
1.5-7 keV, the energy band including the major lines of 
K-shell emissions from Si, S, Ar, Ca, and Fe. We ex- 
tracted the source spectrum from a circular region with 

http: / /heasarc. nasa.gov/docs/suzaku/processing/criteria_xis. html 
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Fig. 1. — Vignetting-corrected XIS image of W49B in the 1.5- 
7 keV band shown on a linear intensity scale. Data from the three 
active XISs are combined. Gray contours indicate every 10% in- 
tensity level relative to the peak surface brightness. The red circle 
and orange annulus indicate the source and background regions, 
respectively. The XIS field of view is shown by the black square. 
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Fig. 2. — Background-subtracted XIS FI spectrum. The energies 
of the prominent emission lines from specific elements are labeled. 

a radius of 2.5 arcmin and the background spectrum from 
an annulus surrounding the source with inner and outer 
radii of 5 and 8 arcmin, respectively. Non X-ray back- 
ground (NXB) spectra constructed with the xisnxbgen 
software were subtracted from both the source and back- 
ground data. After correcting the vignetting and accu- 
mulation area, the source spectrum was made by sub- 
tracting the background spectrum. The spectra of the 
two FI sensors were merged to improve the statistics be- 
cause their response functions were almost identical. Fig- 
ure [2] shows the merged FI spectrum. We can see several 
prominent lines of K-shell emission from He- and H-like 



3.2. APEC Model Fit in the Hard X-ray Band 



As shown in Figure^ the wide-band spectrum includes 
too many emission lines and possibly other complicated 
structures. For clean analyses and discussion, we focus 
on the 5-12 keV energy band, where emission lines and 
other structures mainly originate from Fe. We used only 
the FI data because the NXB count rate of the BI data 
is much larger than the source rate in this energy band. 

We first fitted the spectrum with a VAPEC (CIE 
plasma) model (Smith et al. 2001). The Fe and Ni abun- 
dances (hereafter Z-p^ and ^ni) normalized by the num- 
ber fraction of the solar photosphere (Anders & Grevesse 
1989) were free parameters. Since the Cr and Mn Ka 
emission lines are not included in the VAPEC model, we 
added these lines with Gaussian functions. We fixed the 
interstellar extinction to a hydrogen column density of 
5 X 10^^ cm^^ with the solar elemental abundances, fol- 
lowing Hwang ct al. (2000). To fine tune the energy 
scale, an offset was added as a free parameter and found 
to be -4.7 eV, which was within the allowable range of 
the calibration uncertainties. We added this offset for 
further spectral fitting. Figure [3^ shows the fitting re- 
sult of the one- VAPEC model, with a temperature of 
1.64 keV. The model exhibited significant excess around 
the Fe-Lya line and above 8 keV, and hence was rejected 
with a large x^/dof (degree of freedom) of 1051/138. 

We then tried a two- VAPEC model, assuming equal 
abundances between the two components. Fe-Lya was 
successfully reproduced, but the temperature of one com- 
ponent was unreasonably high ('^70 keV), and a large 
residual above 8 keV remained with an unacceptable 
X^/dof of 562/136. Adding a further plasma or power- 
law component did not improve the fitting any more. In 
every case, the large bump above 8 keV was present. 

To examine whether this bump is real or artificial, we 
checked the light curves of the source and background re- 
gions in the 8.5-10 keV band. They showed almost con- 
stant fluxes, indicating no flare-like event had occurred 
during the observations. In addition, we detected the 
bump in both FI spectra (XIS and XIS 3) and even in 
the BI spectrum. Thus, the bump is a real structure. 

The edge energy of the bump (~9 keV) corresponds 
to the electron binding energy of Fe. This suggests that 
the saw-edged bump is likely due to a radiative recom- 
bination continuum (RRC): X-ray emissions due to the 
free-bound transition of electrons. Hereafter, we call the 
RRC accompanied by the recombination of H-like ions 
into the ground state of He-like ions the Hc-RRC, and 
that of fully ionized ions into the ground state of H-like 
ions the H-RRC. 



3.3. Recombination Continuum and Lines 

Both the line-like excess around Fe-Lya and the bump 
above 8 keV are suggestive of the overionized state be- 
cause there should be a greater fraction of H-like ions 
in the overionized state than in the CIE. This excessive 
number of H-like ions causes strong He-RRC and Fe- 
Lya emissions. However, no current plasma code can 
be applied to the overionized plasma. We, therefore, in- 
troduce the Fe-Lya line and recombination structures 
consisting of the RRC and several emission lines below 
the K-edge energy of the RRC (i?edge), in addition to 
the one- VAPEC model. We consider both He-RRC and 
H-RRC for a consistency check (see Section 4.2). 
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We assume that the RRC is expressed as 



dN 
'dE 



cx exp 



E- E, 



edge 



for E > Ecde, 



(1) 



where kTc is the electron temperature of the relevant 
recombining plasma. This formula gives a good approx- 
imation when the electron temperature is much lower 
than the K-edge energy {kT^ <^ Sedge) (e.g., Smith & 
Brickhouse 2002). In the case of W49B, kT^ is -1.5 keV, 
as we confirm later, while -Eedge of Fe is ~9 keV. Thus, 
we can safely adopt this formula. 

The recombination lines originate from cascade decays 
of electrons that are captured into the excited levels of 
ions by the free-bound transition. We designate the cas- 
cade lines of He-like Fe as Hearec (the principle quantum 
number n = 2— >1) and He/3rec (3— >1). For other lines 
from higher levels {n > 4-^1), we combine them into one 
broad Gaussian function (He7-oO].oc)- The line widths 
and center energies of Hearec and He/?roc are fixed at 
zero and the experimental values for a charge exchange 
process, respectively (Wargelin et al. 2005). On the other 
hand, those of He7 allowed to vary freely. The 

fluxes of all these lines are free parameters. 

The resultant best-fit parameters and models are given 
in Table [1] and Figure [3}d, respectively. With this model, 
X^/doi is greatly improved to 193/128, although this is 
still unacceptable in a purely statistical sense. We can 
see significant data excesses at the lower energy sides of 
the Hea lines of Fe (~ 6.5 keV) and Ni (~ 7.6 keV). 
The former is likely due to the incomplete response func- 
tion. We should note that Fe-Hea statics are superior 
to those for any other objects observed with Suzaku, and 
no significant residual has appeared so far. If we observe 
carefully, we find a similar feature in the Mn-Ka line 
from the onboard calibration source. The latter would 
be partially due to the same reason given above, but is 
mainly due to lack of satellite lines of the Be-like and 
lower ionization states^ of Ni in the VAPEC model. If 
we ignore the 6.4-6.6 keV and 7.6-7.7 keV band to es- 
cape these effects, an acceptable x^/dof of 118/107 is 
obtained. 

4. RESULTS AND DISCUSSION 

We have found, for the first time, the strong He-RRC 
of Fe from W49B. Yamaguchi et al. (2009) recently dis- 
covered the H-RRC of Si and S from a middle-aged 
SNR, IC 443, and hence this is the second discovery of 
a clear RRC from an SNR. We also discovered RRC- 
accompanied recombination lines, which may provide 
good diagnostics for the overionized plasma. We review 
a quantitative verification of our spectral analysis and 
discuss the implications of the results. 

4.1. Contribution of the Recombination Lines 

We discuss the validity of the best-fit fluxes of the RRC 
and recombination lines in Table 1. The recombination 
cross section of the H-like ions into a level of n is approx- 
imately given as shown below (e.g., Nakayama & Masai 
2001). 



1 /3 /cTe 



1 



(2) 



^ As for Fe, such satellite lines are included in the VAPEC model. 
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Fig. 3.— (a) XIS spectrum in the 5-12 keV band. The best-fit 
VAPEC model and additional Ko lines of Cr and Mn with Gaussian 
functions are shown by solid lines. The lower panel shows the 
residual from the best- fit model, (b) Same spectrum as (a), but 
with the radiative recombination continuum of He- like Fe (red), 
recombination lines (orange), and a Lyo line (blue). 

We apply this approximation for the He- like ions, but cti 
must be reduced by half because one electron is already 
at the ground state. In principle, the recombination line 
flux can be estimated by the branching ratio to various 
levels, but these processes are very complicated. We, 
therefore, base our discussion only on a simple picture. 

We compare the predicted capture and observed tran- 
sition rates normalized with the n = 1 value. We can 
estimate the capture rates using Equation 2 as o-2/cti = 
0.62, aa/ui = 0.25, and (cr4 -has -h ...)/ai = 0.34. On the 
other hand, the observed flux rates are Hearec/He-RRC 
= 0.47 (±0.10), He/?rec /He-RRC = 0.053 (±0.023), and 
He7-oOrec/He-RRC = 0.19 (±0.03).^ By taking the frac- 
tions between the observed and predicted rates, we ob- 
tain 76 (±16)%, 21 (±9)%, and 56 (±9)% for the Hearec, 
He/3rec, and He7-c>Oi-oc lines, respectively. 

These fractions for the Hc/Srec and He7-oOrec lines may 
be conceivable, because one electron at n = 1 suppresses 
direct transitions from excited levels {n > 3^1) for the 
He-like ions. The fraction of Hearec is slightly smaller 
than expected because there should be a significant con- 
tribution of the cascade decay electrons from higher lev- 

^ Throughout this paper, all the errors in the parentheses are at 
90% confidence level. 
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TABLE 1 
Best-fit spectral parameters 



CIE (VAPEC) 


feTc [kcVJ 


Z-pc [solarj 


Zni ^solar] 


Normalization^ 




1.52 (1.50-1.53) 


4.44 (4.36-4.53)'" 


10.9 (9.14-12.7)"" 


1.61 (1.60-1.62) 


Line 




Center energy fkcVl 


Line widtli [kcV] 


Normalization'^ 




Cr Ka 


5.655 (5.646-5.663) 


(fixed) 


2.85 (2.59-3.11) 




Mil Ka 


6.162 (6.142-6.183) 


(fixed) 


0.957 (0.750-1.16) 




Fe Hearec'' 


6.666 (fixed**) 


(fixed) 


5.30 (4.20-6.39) 




Fe Lya 


6.961 (6.946-6.971) 


(fixed) 


1.46 (1.27-1.65) 




Fe He/Jrec"^ 


7.880 (fixed«) 


(fixed) 


0.593 (0.338-0.847) 




Fe He7-oOrec'^ 


8.538 (8.492-8.585) 


0.200 (0.163-0.246) 


2.12 (1.81-2.43) 


RRC 




Edge energy [kcV] 


fcTc [kcV] 


Normalization*" 




Fe He-RRC 


8.830 (fixed) 


1.43 (1.30-1.59)' 


11.2 (10.4-11.9) 




Fe H-RRC 


9.194 (fixed) 


1.43 (1.30-1. 59)f 


0.257 (0.00-1.04) 



Note. — The uncertainties in the parentheses are the 90% confidence range. "The unit is 10~^^ / Uenu (iV/(47rD^) (cm~®), where 
Tie, nil, V, and D arc the electron and hydrogen densities (cm~^), emitting volume (cm^), and distance to the source (cm), respectively. 
''These values should be modified in the case of the overionized plasma. See Section 4.3 for details. '"The unit is 10^" cm^'^ s^^. ''Lines 
emitted by a recombination process. '^Fixed at the experimental values (Wargelin et al. 2005). ■'^We assumed the same fcTg. 



els (n > 3^2^1). The real Heaicc flux may be some- 
what larger due to the uncertainty of the response func- 
tion, but the contribution of the Hca^cc flux relative to 
the total Hea is only <10% and does not affect results 
and following discussion. 

4.2. Electron and Ionization Temperatures 

The electron temperatures determined by the 
bremsstrahlung continuum shape and the RRC slope 
are 1.52 (1.50-1.53) keV and 1.43 (1.30-1.59) keV, 
respectively. These consistent results indicate a common 
origin of these emissions. 

The ionization temperatures directly reflect the ion 
fractions and hence can be determined as shown be- 
low. From the best- fit model in Table 1, the flux ra- 
tios of Lya/Hea, He-RRC/Hea, and H-RRC/He-RRC 
are given as 0.016 (0.014-0.018), 0.12 (0.11-0.13), and 
0.023 (<0.10). In Figure 4, we compare these values with 
the modeled emissivity ratios derived from the radiation 
code of Masai (1994) for a plasma of kTc = 1.5 keV. 
We obtain kT^ = 2.58 (2.46-2.68) keV, 2.68 (2.63-2.73) 
keV, and 2.55 (<3.65) keV, respectively, for the above 
ratios. We also compare the Lya/Hca ratio with that 
derived from the APEC code (Smith et al. 2001), al- 
though this code is valid only for a CIE state. We obtain 
kT^ = 2.46 (2.39-2.54) keV, which is within the margin of 
error of the above results. The ionization temperatures 
(~2.7 keV) are significantly higher than the electron tem- 
peratures (~1.5 keV), indicating that the plasma is in a 
highly overionized state. 

The first hint of overionized plasma in W49B was found 
by Kawasaki et al. (2005). Although they analyzed dif- 
ferent elements (Ar and Ca) in different energy bands 
(2.75-6.0 keV), and derived kT^ from the Lya/Hea ratio 
using the CIE plasma code, their results (kT^ ^2.5 keV 
and kTc ~1.8 keV) are nearly consistent with ours. Our 
claim is more essential because it is based on clear detec- 
tion of the recombination structures. 

4.3. Iron and Nickel Abundances 

The abundances listed in Table 1 are valid only for the 
CIE state and should be modified in the overionized case. 
Since no plasma code can be applied to the overionized 
plasma currently, we make possible modifications using 



an available APEC code. 

The Hea intensity is proportional to Z x e(kTc, kT^), 
where Z and e{kTc, kT^) are the abundance of the el- 
ement (solar) and the total emissivity for the Hea for 
kTe and kT^ (cm^s~^), respectively. The emissivities of 
the He-, Li-, Be-, and B-like ions for Fe and the He- 
and Li-like ions for Ni are modified by multiplying the 
ion- fraction ratio between kTz = 2.7 keV and 1.5 keV 
(Mazzotta et al. 1998). The total emissivity is given by 
adding those in individual ionization states. Multiplying 
by e(1.5 keV, 1.5 keV)/e(1.5 keV, 2.7 keV), we obtain the 
real abundances in the overionized state as ZFe~4.9 so- 
lar and ZNi^5.2 solar. Both the elements are highly over 
abundant, indicating an ejecta origin of the plasma. 

4.4. Volume Emission Measure 

To check the consistency of the common origin of the 
bremsstrahlung and RRC emissions, we compare the vol- 
ume emission measure (VEM). The VEM is given by 
/ neniidV/{4:TTD'^), where Ue, nn, V, and D are the elec- 
tron and hydrogen densities (cm^'^), emitting volume 
(cm^), and distance to the source (cm), respectively. 

The VEM of the VAPEC component (VEMyAPEc) is 
derived from Table 1 as VEMvapec = 1-61 (1.60 - 1.62) 
xlO^^ cm-^ On the other hand, the VEM of the RRC 
plasma (VEMyirc) is calculated from 

i^RRC = ai(fcTe) X — X KH-like(A;T^) X VEMkrc, (3) 
riH 

where Frrc, ai(fcTe), npe, and Kn-iikeikTz) are the 
flux of the He-RRC (cm~^s~^), the K-shell recombina- 
tion rate coefficient for kTc (cm'^s^^), the number den- 
sity of Fe (cm^'^), and the ion fraction of H-like Fe for 
kTz, respectively. According to Badnell (2006), the to- 
tal radiative recombination rate of He-like Fe at kT^ = 
1.5 keV is given as ~3.9xl0~^^ cm^s~^. The recombina- 
tion rate into the ground state is given using Equation 2 
as cri/{ai + (72 + •••) ^0.45. The value of npc/nn is cal- 
culated using Z-pc in Section 4.3 and the number density 
of the solar photosphere (Anders & Grevesse 1989). Us- 
ing the observed He-RRC flux, we obtain^ VEMrrc = 

Here, we consider that the kTz error is 2.4-2.8 keV. The large 
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Fig. 4. — Predicted emissivity ratios of Lyo/Hea (black), He- 
RRC/Heo (red), and H-RRC/He-RRC (green) of Fe as a function 
of the ionization temperature (kTz) for an electron temperature of 
1.5 keV (Masai 1994). The horizontal dashed lines represent 90% 
errors of the observed values. 

0.81 (0.58 - 1.60) X 10^3 The two independent es- 

timations of VEM give consistent results, supporting the 
same origin of the overionized plasma in W49B. 

The origin of the overionized plasma in SNRs is an 
open question, and beyond the scope of this paper. We 

error of VEMjijic is due to this effect. 



simply note the possibility of the cooling of electrons 
via thermal conduction (Kawasaki et al. 2005) or a 
more drastic cooling caused when the blast wave breaks 
out of some ambient matter into the rarefied interstellar 
medium (Yamaguchi et al. 2009). If the latter is the 
case, a massive progenitor that had blown a stellar wind 
is likely to be favored. 
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